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Abstract 

QCD as a gauge non-Abelian theory imposes severe constraints on the 
structure of the baryon wave function. We point out that, contrary to 
a widely accepted belief, the traces of baryon number in a high-energy 
process can reside in a non-perturbative configuration of gluon fields, rather 
than in the valence quarks. We argue that this conjecture can be tested 
experimentally, since it can lead to substantial baryon asymmetry in the 
central rapidity region of ultra-relativistic nucleus-nucleus collisions. 
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In QCD, quarks carry colour, flavour, electric charge and isospin. It seems 
only natural to assume that they also trace baryon number. However, this latter 
assumption is not dictated by the structure of QCD, and therefore does not need 
to be true. Indeed, the assignment of the baryon number B — 1/3 to quarks is 
based merely on the naive quark model classification. But any physical hadron 
state in QCD should be represented by a state vector which is gauge-invariant - 
the constraint which is ignored in most of the naive quark model formulations. 
This constraint turns out to be very severe; in fact, there is only one way to 
construct a gauge-invariant state vector of a baryon from quarks and gluons [|IJ 
(note however that there is a large amount of freedom in choosing the paths 
connecting x to Xj): 



B = e ijk 



Pexp (ig J A^dx^ q(x\) Pexp (ig J A^dx^ q(x2 

x Pexp (ig f A^dx^) q(x 3 ) . (1) 

The "string operators" in ([[J) acting on the quark field q(x n ) make it transform 
as a quark field at point x instead of at x n . The e tensor then constructs a local 
colour singlet and gauge invariant state out of three quark fields (see Fig. la). 
The B in eq. ([1]) is a set of gauge invariant operators representing a baryon 
in QCD. With properly optimised parameters it is used extensively in the first 
principle computations with lattice Monte Carlo attempting to determine the 
nucleon mass. The purpose of this work is to study its phenomenological impact 
on baryon number production in the central region of nucleus-nucleus collisions. 

It is evident from the structure of (|l|) that the trace of baryon number should 
be associated not with the valence quarks, but with a non-perturbative config- 
uration of gluon fields located at the point x - the "string junction" |l|]. This 
can be nicely illustrated in the string picture: let us pull all of the quarks away 
from the string junction, which we keep fixed at point x. This will lead to qq 
pair production and string break-up, but the baryon will always restore itself 
around the string junction. The quark composition of this resulting baryon will 
in general differ from the composition of the initial baryon. It is important to 
note that the assignment of the trace of baryon number to gluons is not only a 
feature of a particular kind of the string model, but is a consequence of the local 
gauge invariance principle applied to baryons. 

Surprisingly, at first glance this non-trivial structure of baryons in QCD does 
not seem to reveal itself in hadronic interactions in any appreciable way. This 
is because to probe the internal structure of baryons, we need hard interactions, 
which usually act on single quarks. The baryon number therefore can always be 
associated with the remaining diquark, without specifying its internal structure. 
To really understand what traces the baryon number, one needs therefore to 
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study processes in which the flow of baryon number can be separated from the 
flow of valence quarks. 

In this Letter we suggest that studies of baryon production in the central ra- 
pidity region of ultra-relativistic pp and especially AA collisions provide a crucial 
possibility to test the baryon structure. Our findings may prove important for 
understanding the properties of dense QCD matter produced in AA collisions at 
high energies. In what follows, we shall first formulate our ideas qualitatively, 
and then give some quantitative estimates based on the topological expansion of 
QCD m, H and Regge phenomenology. 

Let us consider first an ultra-relativistic pp collision in its centre-of-mass 
frame, which coincides with the lab frame in collider experiments. At sufficiently 
high energies, the valence quark distributions will be Lorentz-contracted to thin 
pancakes with the thickness of 

where P is the cm. momentum in the collision, and xy ~ 1/3 is a typical 
fraction of the proton's momentum carried by a valence quark. The typical time 
needed for the interaction of valence quarks from different protons with each 
other during the collision is given by the characteristic interquark distance in the 
impact parameter plane, ij ni = const ~ 0(1 fm). However, the time available 
for this interaction in the collision is only t co u ~ zy ~ (xyP)^ 1 . It is therefore 
clear that at sufficiently high energies, when t co u << t int , the valence quarks of 
the colliding protons do not have time to interact during the collision and go 
through each other, populating the fragmentation regions. In the conventional 
picture, the baryon number follows the valence quarks. 

At first glance, the argument looks correct, and is well supported experimen- 
tally - the leading effect for baryons in high energy pp collisions is well established. 
However, the structure of the gauge-invariant baryon wave function (1) suggests 
that this scenario may not be entirely consistent. As we have stressed above, 
in QCD the trace of the baryon number has to be associated with the non- 
perturbative configuration of the gauge field. This "string junction" contains an 
infinite number of gluons which, therefore, by virtue of momentum conservation, 
should carry on the average an infinitely small fraction x$ « xy of the proton's 
momentum. We therefore expect that the "string junction" configuration may 
not be Lorentz-contracted to a thin pancake even at asymptotically high energies 
(see Fig. 2), since 

ZS - » Zy. (3) 

In this case the string junction will always have enough time to interact, and we 
may expect to find stopped baryons in the central rapidity region even in a high- 
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energy collision (of course there will also be baryon-antibaryon pair production 
in addition). This argument leads to a peculiar picture of a high-energy pp 
collision: in some events, one or both of the string junctions are stopped in the 
central rapidity region, whereas the valence quarks are stripped-off and produce 
three-jet events in the fragmentation regions. Immediately after collision, the 
central region is then filled by a gluon sea containing one or two twists, which 
will later on be dressed up by sea quarks and will form baryon(s). Note that 
the quark composition of the produced baryons will in general differ from the 
composition of colliding protons. 

Why then is the leading baryon effect a gross feature of high-energy pp col- 
lisions? The reason may be the following. The string junction, connected to all 
three of the valence quarks, is confined inside the baryon, whereas pp collisions 
become on the average more and more peripheral at high energies. Therefore, 
in a typical high-energy collision, the string junctions of the colliding baryons 
pass far away from each other in the impact parameter plane and do not inter- 
act. One can however select only central events, triggering on high multiplicity 
of the produced hadrons. In this case, we expect that the string junctions will 
interact and may be stopped in the central rapidity region. This should lead to 
the baryon asymmetry in the central rapidity region: even at very high energies, 
there should be more baryons than antibaryons there. 

Fortunately, the data needed to test this conjecture already exist: the exper- 
imental study of baryon and antibaryon production with trigger on associated 
hadron multiplicity has been already performed at ISR, at the highest energy 
ever available in pp collisions || . This study has revealed that in the central ra- 
pidity region, the multiplicities associated with a proton are higher than with an 
antiproton by ~ 10%. It was also found that the number of baryons in the central 
rapidity region substantially exceeds the number of antibaryons ||. These two 
observations combined indicate the existence of an appreciable baryon stopping 
in central pp collisions even at very high energies || . 

Where else do we encounter central baryon-baryon collisions? In a high energy 
nucleus-nucleus collision, the baryons in each of the colliding nuclei are densely 
packed in the impact parameter plane, with an average inter-baryon distance 

r ~ p-^A- 1 '^ (4) 

where p is the nuclear density, and A is the atomic number. The impact parameter 
b in an individual baryon-baryon interaction in the nucleus-nucleus collision is 
therefore effectively cut off by the packing parameter: b < r. In the case of a lead 
nucleus, for example, r appears to be very small: r ~ 0.4 fm, and a central lead- 
lead collision should therefore be accompanied by a large number of interactions 
among the string junctions. This may lead to substantial baryon stopping even 
at RHIC and LHC energies. 
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We shall now proceed to more quantitative considerations. In the topological 
expansion scheme the separation of the baryon number flow from the flow 
of valence quarks in baryon-(anti)baryon interaction can be represented through 
a t-channel exchange of the quarkless junction-antijunction state with the wave 
function given by 



K = e ljk /i' k ' 



( f x2 m * r ( f x2 

Pexp ( ig / A^dx^ \ Pexp ( ig \ A^dx^ 



Pexp (ig I 



A^dx^ 



(5) 



The structure of the wave function (|) is illustrated in Fig. lb - it is a quarkless 
closed string configuration composed from a junction and an antijunction. In 
the topological expansion scheme, the states ([5]) lie on a Regge trajectory; its 
intercept can be related to the baryon and reggeon intercepts |L|: 

a J (0) ~ 2a B (0) - 1 + 3(1 - a R (0)) ~ ^ (6) 

where the baryon intercept «s(0) has been set equal to 0, and the reggeon inter- 
cept or(0) equal to 1/2. 

The Mq exchange should dominate the proton-antiproton annihilation at high 
energies [p]]. Indeed, annihilation requires the baryon number transfer in the 
t— channel. Conventionally, this corresponds to the baryon exchange with the 
intercept ob(0) — 0. Since in Regge theory the energy dependence of the cross 
section is given by s™'^ -1 , and «q(0) > c*b(0), the M(f exchange should give 
the dominant contribution at high energies (see Fig. 3), leading to the following 
energy dependence of the annihilation cross section: 

instead of the s _1 dependence implied by conventional baryon exchange (sq ~ 
1 GeV 2 is the usual parameter of Regge theory). Up to the highest energies where 
the annihilation can still be experimentally distinguished from other inelastic 
processes, the energy dependence (0) is confirmed by the data. Moreover, the 
entire difference between the total pp and pp cross sections can be attributed to 
annihilation (see for a recent review). 

Let us now turn to the consideration of baryon stopping in pp collisions. The 
relevant diagrams are shown in Fig. 4; we consider the simultaneous stopping of 
the two string junctions in the central rapidity region, accompanied by three-jet 
events in the fragmentation regions (Fig. 4a), and the stopping of the junction 
of one proton in the soft parton field of the other, accompanied by one three-jet 
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event (Figs. 4b, c). To calculate the cross sections, it is convenient to evaluate the 
discontinuity of the corresponding three-particle elastic scattering process || [7| 
(see Fig. 5). Introducing the four-momentum ps of the produced baryon (or the 
total momentum of the pair of baryons in the diagram of Fig. 5a) and the four 
momenta of the colliding protons pi, P2 we have the following expressions for the 
invariant energy in the proton-baryon systems: 

sx ~ y/s m t e~ y , 

s 2 ~ \fs m t e y * , (8) 

where s = {P1+P2) 2 is the c.m.s. energy squared of the pp collision, y* is the c.m.s. 
rapidity of the produced baryon(s), and m t is its transverse mass m 2 = m 2 B + p 2 B - 
The product of the invariants (|8]) satisfies the relation 

S1S2 — m 2 s. (9) 

Let us denote the coupling of the Mq reggeon and pomeron to the proton by 
G^ and Gp respectively, and introduce scalar functions fB M ( m t) an d fB P ( m t) 
describing the "two baryons - Mq — Mq" and "one baryon - Mq - Pomeron" 
vertices. The standard calculation |6], [7| then allows us to calculate the cross 
sections; for the diagram of Fig. 5a we get 

E B d -^- = te[G? (0)] 2 f r(m 2 t ) • (10) 

Analogous calculation for the sum of diagrams in Figs. 5b, c gives 

= B*G? (0)GJ(0) /fK) 

x (exp[y*(a P (0) - a J (0))} + exp[-y*(a P (0) - a J (0))}) . (11) 

Using the value (f| of the Mq intercept, we find that the double baryon produc- 
tion cross section (10) has ~ s -1 / 2 energy dependence and (within the central 



rapidity region, where our considerations apply) does not depend on rapidity. 
The cross section (|ll|) of single baryon stopping also decreases with energy, but 
much more slowly. Writing down the Pomeron intercept as 

a P (0) = l+A, (12) 

one gets the energy dependence of the cross section ( |TTD in the form ~ s -i/4+a/2 
At very high energies the process of single baryon stopping will therefore be more 
important. Note however that the ISR data show a large value of the correlation 
between the probabilities of the stopping of the beam baryons, indicating that 
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the process of double baryon stopping may still dominate even at rather high 



energies. The rapidity dependence of (|TT|) does not show a "central plateau", 



indicating instead a "central valley" structure. This is in accord with the ISR 
data; moreover the particular rapidity dependence of (pTT|) reproduces the data 
H reasonably well. 

The functional dependence of the single baryon stopping cross section sim- 



ilar to (|il|) has been advocated before || in a different approach; the authors 
considered a specific mechanism of perturbative destruction of the fast diquark, 
accompanied by the baryon number flow over a large rapidity gap. In the frame- 
work of their approach, the authors of ref.[[J have also performed a calculation 
of the cross section, based on the combination of perturbative technique and 
constituent quark model. 

Unfortunately, since we believe that the dynamics of the stopping process is 
genuinely non-perturbative, so far we have not been able to find a reliable way of 
computing the constant G^(0) entering the expression ([TT|) and can only extract 
it from the existing data. However once it is done, we can perform parameter-free 
extrapolation to higher energies. 



As we have already stressed above, the formulae ( |l~0| , |i~l"|) refer to the net baryon 
number, i.e. they refer to the difference between the number of produced baryons 
and antibaryons. The process of baryon-antibaryon pair production will therefore 
represent an important background to baryon stopping; the data |§ show that 
at ISR energies the probability of baryon stopping is about three times smaller 
than the probability of baryon pair production. At high energies, the dominant 
contribution to the BB pair production will be given by the interaction of two 
Pomerons, with the cross section 

rJ 3 rr( BB ) ( R m \ 2q H0)-2 

e b-^7- = ^Hm 2 /IbK) , (13) 

representing an energy-independent fraction of the total cross section: 

/ a \ «p(0)- 1 

a BB ~ a toi ~ . (14) 

Since the topological structure of the Pomeron is that of a cylinder [[l ] 

P = Tr 



Pexp yig j> A^dx^ 



(15) 



at high energies the associated multiplicities in the processes of single (n^) and 
double [rS 2 ^) baryon stopping will be higher than in the average inelastic event, 
described by the cut of the Pomeron exchange diagram (see Fig. 6): 

n (l) ~ 1 n inel. n (2) _ _ tfnel ^ (jg) 
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as we already discussed above at the qualitative level. Also, since the baryon 
stopping and baryon pair production arise in our scheme from different kinds of 
t-channel exchange (M(j and the Pomeron, respectively), we expect a small corre- 
lation between the events with the pair production and stopping. Experimentally, 
it was found to be 0.16 ± 0.22 §. 

Before we turn to the discussion of nucleus-nucleus collisions, it is useful to 
recall the geometrical picture of high energy scattering in the impact parameter 
plane. In the impact parameter representation, the growth of the total cross 
section at high energies as described by one-Pomeron exchange can be attributed 
to the increase of the effective radius of the interaction, according to 



R mt ~^2Rl + a> P ln(s/s ), (17) 

where R p is a constant and a' P is the slope of the Pomeron trajectory. (The 
Froissart bound allows even faster growth of the interaction radius with energy: 
R%nt ~ ln(s/so).) The central region of the disk becomes completely black at 
high energy if op(0) > 1; a further growth of the interaction strength in the 
centre of the disk is prevented by the unitarity constraint imposed on the partial 
amplitudes. 

The colliding nuclei, the transverse plane nucleon distributions in which are 
characterized by the packing parameter (|J) will therefore see each other as uni- 
form black disks. This means that in a central nucleus-nucleus collision the cross 
section of the inelastic nucleon-nucleon collision will not further increase with 
energy when Ri nt » r since the soft peripheral interactions building up the 
Pomeron will be effectively screened out. On the other hand, the processes of 
baryon stopping are central in the impact parameter plane, and therefore may 
not be screened in the case of nuclear collisions. A very slow decrease of the cross 
section fllTD with energy implies then that even at LHC energies the nuclear 
stopping may still be present, as we shall now discuss. 



The ISR data |10|, [TT[] show that at y/s = 53 GeV the cross sections of proton 



and antiproton production at y* = and pt = 0.6 GeV '/c are 

= 0) = 0.700±0.162 mb GeV~ 2 ; ^(y* = 0) = 0.430±0.033 mb GeV~ 2 . 
d A p d A p 

(18) 

Since the mechanism of baryon- antibaryon pair production (see (0)) obviously 
leads to an equal number of protons and antiprotons, the data fll8Q imply that 
the following fraction of the protons in the central rapidity region is produced by 
stopping: 

UW~s = 53 GeV) = °^£= ~ 25%. (19) 
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We can now estimate the ratio R of multiplicities associated with proton (n p ) and 
antiproton (nP) production, using the predictions fll6|) . Since the multiplicity 
associated with protons and antiprotons should be the same in the absence of 
stopping, we get 

R=- = (l-f st )+f st (-), (20) 
where we have omitted the contribution of the double stopping process, which 



is asymptotically suppressed at high energies according to (^0|) , but may still 
be important at ISR energies ||. Assuming that the multiplicity associated 
with antiprotons does not differ substantially from the average multiplicity of an 
inelastic event, nF ~ n mel , in accord with (|H|) and with experimental data ||, 
we obtain from ( p0|) and ([Tj]) an estimate 

#-1.05. (21) 

Even though the value fl2"T|) agrees within experimental errors with the measured 
value of R ~ 1.1 0, we may conclude that the contribution of double baryon 
stopping with higher associated multiplicity ( |IED is possible. A detailed exper- 
imental study of double baryon production in pp collisions would therefore be 
useful to clarify the situation. 

Extrapolation to the pp collisions at the energies of ^/s ~ 6 TeV (correspond- 
ing to the c.m.s. energy per nucleon-nucleon collision in Pb-Pb interactions at 
LHC) according to formulae (|TT| , |13| ) with A = 0.08 [|l^] yields then the following 



stopping fraction: 

f st {^ = 6 TeV) ~ 5%. (22) 

The nucleus-nucleus collisions will be accompanied by much larger stopping, as 
we discussed above, and we expect that the estimate (E21) in this case can only be 



considered as a lower bound on the baryon asymmetry. Therefore we may expect 
substantial excess of baryons over antibaryons in the central rapidity region of 
nucleus- nucleus collisions at LHC. 

It is interesting that already at SPS energy, the otherwise successful phe- 
nomenological approaches based on the topological expansion [13-16], but not 
taking into account the presence of the string junction explicitly, seem to under- 



estimate [17] the pronounced baryon stopping observed experimentally in nucleus- 



nucleus collisions ^M. 

It would be useful to analyse the dynamics of baryon stopping in an approach 
where the topological structure of the baryon is explicit: the Skyrme model. The 
formation of baryon-antibaryon pairs in this approach is treated as the forma- 
tion of topological defects in the quark condensate |19| (the net baryon number 
of the produced pairs is of course equal to zero). The picture proposed above 
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would invoke into consideration also the high-energy scattering of such topolog- 
ical defects. Since the Lorentz boost of such configurations is not trivial, one 
may expect the occurence of the final states in which one or two Skyrmions are 
stopped in the central region, and the part of their pion field is "shaken off" to the 
fragmentation region. Such events would contribute to the baryon asymmetry in 
the central region. We leave the consideration of the baryon stopping dynamics 
in topological models for further studies. 

Amazingly, the so-called "Centauro" and "Chiron" events reported in cosmic- 



ray emulsion experiments [2(| and interpreted in favour of the existence of dense 



quark matter by Bjorken and McLerran |2l| , are characterized by non- vanishing 



baryon number density. As was formulated in ref. ||22| , "the fluid in the central 
region has no net baryon number, so that there would need to be a spontaneous 
generation of net baryon density to make these objects". Indeed, in the scenario 
proposed by Bjorken fl22| , at ultra-relativistic energies the valence quarks of the 



colliding nuclei pass through each other, leaving behind a "little Universe" of zero 
net baryon density, which at the moment of its production contains a gluon sea. 
In our picture, this sea from the very beginning is made stormy by the presence 
of non-perturbative twists - specific configurations of the gluon field, which trace 
non- zero net baryon number. The "little Universe" , just like our big one, may 
therefore generate substantial baryon asymmetry. 

I am very indebted to G. Veneziano for many instructive suggestions and 
encouragement. It is a pleasure to acknowledge stimulating and enlightening 
discussions with J.-P. Blaizot, A. Capella, J. Ellis, K.J. Eskola, M. Gazdzicki, 
K. Kajantie, L. McLerran, J.-Y. Ollitrault and H. Satz, whom I also thank for 
his interest in this study. This work was supported by the German Research 
Ministry (BMBW) under contract 06 BI 721. 
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